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We have measured the specific heat of an S = 1/2 antiferromagnetic alternating Heisenberg chain 
pentafulorophenyl nitronyl nitroxide in magnetic fields up to H > Hc2- This compound has the 
field-induced magnetic ordered (FIMO) phase between Hci and Hc2- Characteristic behaviors are 
observed depending on the magnetic field up to above Hc2 outside of the H-T boundary for the 
FIMO. Temperature and field dependence of the specific heat are qualitatively in good agreement 
with the theoretical calculation on an S = 1/2 two- leg ladder. [Wang et al. Phys. Rev. Lett 84 
5399 (2000)] This agreement suggests that the observed behaviors are related with the low-energy 
excitation in the Tomonaga-Luttinger liquid. 



Since the Haldane conjecture [l|, the quasi-one- 
dimensional (ID) antiferromagnetic spin-gapped sys- 
tems, have been attracting much attention from both ex- 
perimental and theoretical points of view. Spin-gapped 
systems, such as Haldane systems, spin-Peierls systems, 
two-leg ladder systems and alternating chain systems 
have a finite energy gap A between the non-magnetic 
ground state and the lowest excited state because of no- 
table quantum effects. When magnetic field is applied 
to these systems, the lowest branch of the first excited 
state goes down due to the Zeeman effect and intersects 
the ground state at the lower critical field ifci=A/#/iB, 
indicating that the gap closes. For the higher field, the 
magnetization starts to grow and saturates at the up- 
per critical field Hq2 with a spin-polarized state. For 
Hci < H < Hq2: the system has the gapless excita- 
tion on the magnetic ground state, and is described by 
the universarity class of ID quantum systems called the 
Tomonaga-Luttinger liquid (TLL), in which heuristic fea- 
tures, such as the power-law temperature dependence of 
the NMR relaxation rate and the linear-T dependence of 
the low-temperature specific heat are theoretically pre- 
dicted 2]. Experimentally, Chaboussant et al. reported 
the observation of the TLL state in an S = 1/2 two-leg 
ladder Cu 2 (C 5 Hi 2 N 2 )Cl4 (CuHpCl) in the NMR study 
in magnetic fields up to above Hq2 L| • As for the previ- 
ous specific heat measurements on the spin-gapped com- 
pounds, however, most of those results were limited far 
below Hq2 and the probing of the TLL behavior by mea- 
suring the specific heat in the entire field range of the 
gapless phase has not been done [4, 5, 6-8]. 

The compound we focus on here is a genuine or- 
ganic radical crystal pentafulorophenyl nitronyl nitrox- 
ide (F 5 PNN). The magnetism of F 5 PNN comes from 
an unpaired electron localized around NO groups in the 
molecule. The low- temperature susceptibility is well re- 
produced by an alternating chain model with alternating 
exchange interactions along the chain, Ji/ks = —2.8 K, 
J2/J1 = 0-4 j2|. This compound has two- fold advan- 
tages. Firstly, it provides us with a really isotropic quan- 



tum spin system because of the quenching of the spin- 
orbital interactions. In other words, we are free from 
the anisotropy effects which are ordinarily inevitable in 
inorganic compounds. Secondly, this compound has an 
appropriate values of Hci and Hq2 for the experimental 
study. From the low-temperature magnetization process, 
Hci and Hc2 were estimated to be 3 T and 6.5 T, re- 
spectively 9]. It should be noted here that these two 
critical fields are much smaller than those of other spin- 
gapped compounds and can be easily attained with a 
conventional superconducting magnet, which enables us 
to perform the detailed specific heat measurements on 
F5PNN in the three different phases. 

In our previous paper, we have found the field-induced 
magnetic ordering (FIMO) in F5PNN by the specific heat 
measurements in magnetic fields up to 6 T and at around 
200 mK 10]. At H = 5 T, a characteristic hump was 
observed at around 500 mK outside of the FIMO phase 
in the H-T phase diagram. On the other hand, Izumi et 
al reported that the NMR relaxation rate Tf 1 of F 5 PNN 
showed an anomalous power-law behavior below about 4 
K for Hci < H < Hc2 From the experimental point 
of view, however, we have to be careful about the effect 
of the inevitable inter-chain coupling J' on the present 
physics of interest. In F5PNN, J' jk& is estimated to 
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FIG. 1: Magnetic field versus temperature phase diagram 
for F 5 PNN based on data from the previous [T(| and present 
specific heat measurements. 
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be 1 x 10 2 K by the mean-field theory for FIMO in 
the gapless phase, where the three-dimensional ordering 
temperature Tq is given as Tq w (Here we take 

J/A;b as the average value of Ji/fce and ^/ks and Tq — 
200 mK at if = 5 T for an order estimation). The value 
of J'//cb gives an exchange field if ex c± 2J g ^ of about 
10 2 Oe, about 0.5 % for the external field 5 if = 5 T. 
Thus, the hump in the specific heat observed outside of 
the FIMO boundary is intrinsic to the one-dimensionality 
of F5PNN, and should be studied more in detail. 

In this letter, we extend the specific heat measurement 
in magnetic field up to if = 8 T (> if 02) an d in the 
temperature range up to 6 K. In addition to the sharp 
anomaly due to FIMO, a hump or a shoulder in the spe- 
cific heat is observed in the field range above if ci- The 
temperatures of these anomalies are plotted in if -T phase 
diagram of Fig. 1, together with our previous results [HI- 
Although the low-temperature specific heat at the crit- 
ical field if ~ if ci shows linear-T dependence, it does 
not seem to approach zero monotonously as T — > K 
at if = 2.5 T and 7 T, where the system is just outside 
of the gapless phase. These results will be discussed in 
connection with the TLL behaviors which were reported 
in the above NMR study and in the theoretical specific 
heat study of the two- leg ladder system [l^ . 

The sample preparation of F5PNN is detailed in Ref. 
[l^ . Specific heat measurements on the polycrystalline 
sample of F5PNN were performed by the adiabatic heat- 
pulse method using a 3 He- 4 He dilution refrigerator. 

Figure 2 shows some representative data for tempera- 
ture dependence of the total specific heat C in the mag- 
netic fields up to 8 T. First, we give characteristic be- 
haviors of C in some typical fields. At if = T, C 
shows a broad maximum around T = 2 K, which is due 
to the short-range ordering in a ID magnetic system. 
At low temperatures, C decreases exponentially with de- 
creasing temperature, reflecting the energy gap above the 
non-magnetic ground state. As magnetic field increases, 
the broad maximum is suppressed and broadened. In 
the field if = 3 T(~ ifci)? C shows linear-T depen- 
dence. At 5 T, a sharp anomaly due to the FIMO is 
seen at T ^ 200 mK and the hump appears at around 
500 mK with its maximum. These results are consistent 
with our previous work. In magnetic fields higher than 
if C2(— 6.5 T), C shows a shoulder in the temperature 
range of 500 mK< T < 1.5 K and decreases exponen- 
tially with decreasing temperature. In order to see these 
features more clearly, we estimated temperature depen- 
dence of the magnetic specific heat C m by subtracting 
the lattice contribution C\ from C. A broken line in Fig. 
2 shows C\ which is estimated from the data at if = 
T so that the total magnetic entropy for N spins should 
approach Nk^hi(2S J r 1) at higher temperatures. C m will 
be shown in the following figures in the four field regions 
covering the two critical fields if ci and if 02- 
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FIG. 2: Some representative data for temperature depen- 
dence of C of F5PNN in magnetic fields up to 8 T. A broken 
line shows the lattice contribution. 



In Fig. 3(a), we show C m in the field range of if < 
2.5 T. All data exhibit the exponential decay characteris- 
tic of the spin-gapped phase at low temperatures. With 
increasing magnetic field, the exponential behavior grad- 
ually changes to the almost linear-T dependence, indicat- 
ing the suppression of the energy gap. It is noteworthy 
that C m does not seem to approach zero monotonously 
as T ^ K at ii = 2.5 T, which may imply a certain 
structure in C m at lower temperatures. Because 2.5 T is 
lower than ifci, it is hard to consider that this behavior 
is relevant to the FIMO. We will discuss this point later. 

Figure 3(b) shows the results of C m in the field range 
of 3.0 T < if < 6.5 T, which corresponds to the gapless 
phase between if ci and ifc2 estimated by the magnetiza- 
tion study Q • The characteristic behaviors in this phase 
are emergence and shift of the hump and appearance of 
the shoulder depending on the field strength outside of 
the FIMO phase in the if-T phase diagram. At 3 T, 
an up-turn deviating from linear-T dependence is seen 
at lowest temperatures, which can be due to the FIMO. 
The up-turn is enhanced with increasing field up to 5 T, 
and seems to be suppressed in the higher field. At if = 5 
T, it is clear that the hump emerges at around 500 mK 
as indicated by the down- arrow, followed by the sharp 
anomaly due to the FIMO at T ~ 200 mK as mentioned 
before. Further, a shoulder emerges at around T ~ 700 
mK at if = 6.5 T as indicated by the up- arrow. 

For a closer look of the results in this region, the low 
temperature parts of C m for4T<if<6.5T are shown 
in Fig. 3(c). At if = 4.5 T, the hump appears at T ~ 
600 mK. As increasing field, the hump shifts to the lower 
temperatures. Finally, the hump merges with the sharp 
anomaly due to the FIMO at if = 5.5 T and is not ob- 
served at if = 6 T. From the emergence and shift of 
the hump, it is reasonable to consider that the hump is 
masked by the broad maximum due to the short-range 
ordering in the field range of 3 T< if < 4 T. While the 
hump already has not been seen at all at 6 T, the shoulder 
emerges at 6.5 T at around 700 mK above the temper- 
ature range where the sharp anomaly due to FIMO is 
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FIG. 3: Temperature dependence of C m in magnetic fields, (a) T< H <2.5 T (spin-gapped phase) (b) and (c) 3 T< H <6.5 
T (gapless phase) (d) 7 T< H <8 T (spin-polarized phase). A solid line in (b) shows linear-T dependence of C m - Down-arrows 
in (b), (c) indicate a hump in C m . Up- arrows in (b), (c) and down-arrows in (d) indicate a shoulder C m . Insets show theoretical 
calculations of the strongly coupled two- leg ladder system [T^. 



seen. The origin of the behaviors will be discussed later. 

Figure 3(d) shows temperature dependence of C m for 
H>7 T, where the system is in the spin-polarized phase. 
We can see the exponential decay at low temperatures 
as in the case for H<2.b T. In this phase, the spins 
are forced to align parallel by a strong magnetic field. 
That is, the decay indicates an energy gap induced by 
the strong magnetic field. The shoulder in C m is also 
seen in the temperature range of 600 mK< T < 1 K as 
in the data at 6.5 T in Figs. 3(b) and (c). The magni- 
tude of this shoulder at 7 T is larger than that at 6.5 T 
and locates at higher temperatures. This shoulder fur- 
ther develops with increasing magnetic field. It is noted 
that C m does not seem to approach zero monotonously 
down to K at H = 7 T implying a certain structure in 
C m at lower temperatures as in the case at H = 2.5 T. 

To discuss the present results, we first refer to a theo- 
retical study by Tachiki and Yamada about the specific 
heat for weakly coupled 5=1/2 dimer systems, which in- 
clude the alternating chain system and the two-leg ladder 
system in certain parameter regions, in magnetic fields 
(l4j . In the study they obtain an effective Hamiltonian in 
magnetic fields around Hqi by neglecting the upper two 
levels in the energy spectrum of each dimer and replacing 
the lower two levels, corresponding to the singlet ground 
state and the lowest branch of the triplet, by the effec- 
tive spin S = 1/2. This replacement allows us to describe 
the systems by the Hamiltonian of a ID S = 1/2 XXZ 
model with easy-plane anisotropy in the effective field. It 
is noted that physical properties of the alternating chain 
system and the two-leg ladder system in magnetic fields 
are expected to be qualitatively same from this mapping 
of the Hamiltonian. Indeed, the specific heat of F5PNN 
is similar with that of an S = 1/2 two- leg ladder system 
CuHpCl, above the critical temperature of the FIMO [6- 
8] . A hump in C m is seen in the gapless phase of both 
compounds. Whereas, the NMR relaxation rate T^ 1 of 
CuHpCl shows an anomalous increase with the power-law 
temperature dependence as T^O K, which is character- 



istic of the TLL, in the magnetic field and temperature 
range where the hump in C m was observed |2j. In the 
case of F5PNN, the power-law behavior in T^ 1 is also 
observed below about 4 K in the gapless phase, as men- 
tioned before 11]. This temperature-field range includes 
that where the hump in C m is seen. 

Then, we refer to a transfer-matrix renormalization 
group study by Wang et at. for thermodynamic proper- 
ties of an S = 1/2 two- leg ladder model in magnetic field 
with parameters Jm = 1 and J±/J\\ — 5.28, where Jm and 
J± are the intra- and inter-chain coupling, respectively 
[13. In this model, the two critical fields are estimated 
as H C i = 4.3823 and H C 2 = 7.2800. Our results are 
qualitatively in good agreement with their calculations 
on the whole phases. Their calculations for H < Hqi, 
Hqi < H < Hq2, and H > Hq2 are shown in the insets 
of Figs. 3(a), (b) and (c), (d) respectively. In the fol- 
lowing, we will analyze our results by comparing them to 
their calculations in each phase. 

In the spin-gapped phase below Hqi (Fig. 3(a) inset), 
their calculations show a single peak structure due to the 
short-range ordering and suppression of the energy gap 
by applying magnetic field. It should be noted that their 
results show a cusp just below Hqi at very low tempera- 
tures. Turning now to the present experimental results, 
the data at 2.5 T, where is just below Hqi , do not seem to 
approach zero monotonously down to K as mentioned 
before. From the good qualitative agreement between 
both results, it is considered that the above behavior at 
2.5 T is ascribed to the cusp in their calculations. There- 
fore, we will clearly observe the similar cusp if we can 
make the measurement at lower temperatures possible. 
Wang et al. suggest that this cusp appears at the vicin- 
ity of boundary between the spin-gapped and TLL phase 
in the H-T phase diagram obtained from the calculations 
of the specific heat and the magnetization. 

In the field range of Hqi < H < Hq2 (the gapless 
phase : Figs. 3(b) and (c) insets), it is the most signif- 
icant feature in their calculations that a hump emerges 
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and then shifts to low temperatures as increasing mag- 
netic field and a shoulder emerges at H ~ Hq2- This 
feature is also seen in our results as mentioned before. 
We estimated the ratio of the temperature of the hump 
with its maximum to that of the broad maximum due 
to the short-range ordering in the ID magnetic systems 
in both our results and theoretical results. From this 
estimation, the ratio in our results is roughly consistent 
with that in theoretical results. (In both cases, the ra- 
tios are about 10 % in the middle field of the gapless 
field region.) According to Wang et al., the hump comes 
from low-energy excitations in the TLL and the linear- 
T dependence should be observed at low temperatures. 
However, the linear-T dependence in C m is not seen in 
our results because it is masked by the anomaly due to 
FIMO at lower temperatures. As for the origin of the 
shoulder, they do not make a definite remark but sug- 
gest that it appears at the vicinity of the TLL phase. 

In the spin-polarized phase for H > Hq2 (Fig. 3(d) 
inset), the development of the shoulder is clearly seen 
in their calculations as well as our experimental results. 
Moreover, another cusp is seen at low temperatures just 
above Hq2- In the corresponding field of our results 
{H = 7 T) , C m does not seem to approach monotonously 
zero as T^O K. From the good agreement between both 
results, our result at 7 T can correspond to this cusp as in 
the case at 2.5 T. This cusp also locates at the vicinity of 
boundary of the TLL phase in their H-T phase diagram. 

Finally, to evaluate the validity of the above compari- 
son, we will show the field dependence of the maximum 
magnetic specific heat C™ ax and the corresponding tem- 
perature T max of the present results in Fig. 4. As in- 
creasing magnetic field, C™ ax first declines and shows a 
minimum around the field value of (ifci+-Hc2)/2. Then, 
C™ ax begins to increase for the higher fields up to 8 T. 
On the other hand, T max keeps its original value and be- 
gins to increase at H ~ Hqi and reaches a maximum 
around Hq2- These features of field dependence of C™ ax 
and T max well agree with the theoretical results by Wang 
et al. (Fig. 6 of Ref. 12). A theoretical calculation 
is required for the present alternating Heisenberg chain 
system with the experimental features mentioned above. 

In conclusion, we have performed the detailed specific 
heat measurement on the S =1/2 alternating Heisenberg 
chain system F5PNN in magnetic fields up to 8 T. This 
is the first detailed experiment on the specific heat of the 
spin-gapped system up to above Hq2- The characteristic 
behaviors are observed at around the gapless field region 
outside of the FIMO phase, which well agree with the 



theoretical calculation in the TLL regime qualitatively 
[l^ . Thus, this agreement suggests that the observed 
behaviors are related with the low-energy exicitation in 
TLL state. This scenario is also supported by the power- 
law behavior in the NMR relaxation rate T^ 1 of F 5 PNN 
observed below 4 K in the gapless field region [l]| . 

We would like to express our sincere thanks to A. 
Tanaka, K. Izumi, T. Goto, S. Suga, N. Maeshima and 
T. Sakai for valuable discussions. 




FIG. 4: The maximum magnetic specific heat C™ ax and the 
corresponding temperature versus magnetic field H. Er- 
rors at around Hq2 are due to the broader maximum of C ma g • 
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